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Abstract. Expressions for surface brightness distribution and for flux density have been 
theoretically derived in the case of two simple models of a shell supernova remnant. The 
models are: a homogenous optically thin emitting shell with constant emissivity and a 
synchrotron shell source with radial magnetic field. Interactive Data Language (IDL) codes 
for fitting theoretically derived emission profiles assuming these two models to mean profiles 
of shell supernova remnants obtained from radio observations have been written. 

1. MODELS OF EMISSION OF SHELL SUPERNOVA REMNANTS 

1. 1. HOMOGENOUS OPTICALLY THIN EMITTING SHELL WITH CONSTANT EMISSIV- 
ITY 

In this paper, we investigated two models of emission of shell supernova remnants 
(SNRs). If we consider homogenous emitting shell with emissivity e v — const, for 
specific intensity, if the medium is optically thin, we have 

e » dS -\e„(r> 2+ -rU, ^<sin0<f, (1) 

where ds = dr'. Cosine theorems (see Fig. 1) 

(R - A) 2 = d 2 + rf - 2dr[ cos 0, (2) 
R 2 = d 2 + r 2 2 - 2dr' 2 cos 0, (3) 

give us 



r' 1± = d cos 9 ± J (R - A) 2 - d 2 sin 2 0, (4) 



'2± 

Finally, we have 



rU = dcos9± VR 2 - d 2 sin 2 9. (5) 



C v (Vsin 2 6» 2 - sin 2 9 - ^sin 2 9 l - sin 2 6>) , < 9 < 9 X 
CWsin 2 2 - sin 2 0, 9 l < 9 < 9 2 



(6) 



1 



S. OPSENICA and B. ARBUTINA 




Figure 1: Radiation from an optically thin homogenous shell with thickness A and 
radius R, at the distance d from the observer. r[ and r' 2 are the distances from 
intersection of the line of sight for a given 9 with the inner and outer radius of the 
shell, respectively. 

where 9\ — arcsin R ~ d A , 62 = arcsin ^ and C v — 2e v d. From the last equation one 
can see that J° = 2e v RS and /™ ax = 2e u R v / S(2 - S), where S = A/R. Brightness 
distribution i.e. specific intensity in units 2e v d given by equation (6), for 5 = A/R = 
0.1 and R/d = 0.01, can be seen in Fig. 2. 
For the total flux density we have 

S v = / Iv cos ^ sm ^d^dy) 

Jo Jo 

= 47r ^^ 1 (v 7 ^ 2 -d 2 sin 2 6»- \J{R— A) 2 - d 2 sin 2 fl) cos 61 sin 6d0 
+ 4tts v [ VR 2 - d 2 sin 2 9 cos 9 sin <9d6». (7) 
After integration we obtain the expected result 

where the shell volume is V — / = 1 — (1 — (5) 3 is the volume filling factor, 

£ v = A-KEu is total volume emissivity [e v is emissivity per unit solid angle) and L v is 
luminosity. 

1.2. SYNCHROTRON SHELL SOURCE WITH RADIAL MAGNETIC FIELD 

If we have a synchrotron shell source with radial magnetic field, emission coefficient 
is e v oc {Bsm6') a+1 is- a i.e. 

e v = i„{sm6') a+1 - (9) 
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Figure 2: Brightness distribution for an optically thin homogenous shell-like source 
with S = A/R = 0.1 and R/d = 0.01. 

Sine theorem (see Fig. 1) gives us 

7 = ^1. Q = 9' - 9, (10) 
a sin 0' 

r = d(cos9-sm9 cot 9') (11) 

ds = dr'=dsm9 =— (12) 

sin 9' 



i.e. 



and 



Intensity is then 



I v = J e v ds = e v dsm9 j {sm9') a - 1 &9 1 



i.e. 



(13) 



] 2C V sin 9 (1 - < 9 < 9, 

' C v sin 9 /; 2 2 ; (1 - M 2 ) ^ d/i, 0i < < 2 , 



sin Si j2 

Rather than direct integration we will find flux density through S v — — . 



where /x = cos9', /Ui,2± = T "^'"sil^ 2 S '" ~ an d = 

Ra 
where 

£ v = e„du' = / e v (sm9') a+1 sm9'd9'dtp = 2ni v / (sin6»') a+2 d6l' (15) 

J47r JO JO Jo 
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Figure 3: Radial profile with S = A/R = 0.1 and R/d — 0.01 convolved with a 
Bessel function with HPBW=5 x 10~ 4 rad. 

1. e. £ v = 27ry / 7r^^r|g,y and the shell volume is as before V = / = 1 — (1 — <5) 3 . 

2. IDL CODES 

2. 1. DIRECT PROBLEM: SIMULATION OF OBSERVATIONS OF SHELL SUPERNOVA 
REMNANTS WITH RADIO TELESCOPE 

If we observe a shell SNR with a radio telescope, picture we get is a convolution of 
real intensity of radiation of the SNR and power pattern of the telescope, so we get 
"convolved" intensity (Fig 3). When simulating this convolution numerically, one 
must choose an expression for real intensity and an expression for power pattern of a 
radio telescope. In our case, expressions for intensities are (6) and (14), according to 
the models. For power pattern P n (9), two possible cases have been chosen: Gaussian 
approximation and approximation with Bessel function of the first kind. Each pat- 
tern has defined half power beam width (HPBW). Usually one takes HPBW=1.02-^ 
from Bessel function approximation (D is diameter of radio telescope and A is wave- 
length, see Rohlfs and Wilson 1996, Urosevic and Milogradov- Turin 2007). Because 
of technical limitations, we must consider that power pattern takes zero value for 
angles greater than some critical angle 9 C . In the case of Gaussian approximation of 
power pattern, 9 C of 5 sigmas (cr = HPBW/(2\/ln 2)), while for the approximation 
with Bessel function 9 C of 8 HPBW has been chosen. 

Expression that is used for numerical simulation of convolution of intensity of 
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radio emission from a SNR and power pattern of a radio telescope is: 

/C o„v , ) = I /intersection ^ W P " ( g ' ) Sm 0d0d( f 

2tt J° c P n (9) sin 9d9 

Angle 9 is related to other angular parameters through following relation of spherical 
trigonometry: 

cos 9 — cos 9 cos 6*o + sin 9 sin # cos ip. (17) 

Region of double integration in numerator of the expression (16) is the intersection 
between regions where two convolving functions I v (9) and P n (9) are defined. That 
double integration is performed by the IDL function INT_2D. Integration in denomi- 
nator of the expression (16) is performed by the IDL function QROMB. Integrations in 
the expression (14) are performed by the IDL functions QROMB and I NT_TAB UL ATED , as 
well as by "handwritten" function that calculates definite integrals using rectangular 
method. 

In the case of first model (with constant emissivity), user of the program enters 
parameters C v , 9\, 9 2 , as well as parameter of antenna HPBW, and the program 
performs a convolution. In the case of second model (with radial magnetic field) , user 
also enters an additional parameter of object - spectral index a. 

2. 2. INDIRECT PROBLEM: FITTING MODEL TO OBSERVED PROFILE OF A SUPER- 
NOVA REMNANT 

Indirect problem is the following: user enters observed radial emission profile of a shell 
SNR in the form of a table, as well as the parameter HPBW, and the program should 
find the best values for parameters C„, 9\ and 9i in the case of first model, or C v , 
9\, 9i and a in the case of second model, by fitting the chosen model to the entered 
data. This is performed by the iterative IDL procedure CURVEFIT. To perform this 
procedure, user has to estimate initial values for the parameters. That can be done 
using observed radial profile and the expression (6) or (14). Initial value for spectral 
index a can be taken to be 0.5. This parameter is, however, better to be held fixed 
(assuming that it is known from spectra). In addition to finding the best values of 
the parameters, the program also calculates their errors (i.e. standard deviations). 
Finally, the program calculates the flux density of an SNR. 

The program has been tested with artificially generated data. Results of the 
application of the program to the real data will be given elsewhere. 
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